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We study experimentally and theoretically the effect of Eu doping and partial oxygen iso- 
tope substitution on the transport and magnetic characteristics and spin-state transitions in 
(Pri_yEu y )o.7Cao.3Co03 cobaltites. The Eu doping level y is chosen in the range of the phase 
diagram near the crossover between the ferromagnetic and spin state transitions (0.10 < y < 0.20). 
We prepared a series of samples with different degrees of enrichment by the heavy oxygen isotope 
ls O, namely, containing 90%, 67%, 43 %, 17 %, and % of ls O. Based on the measurements of 
ac magnetic susceptibility x(^) an d electrical resistivity p(T), we analyze the evolution the sample 
properties with the change of Eu and 18 O content, ft is demonstrated that the effect of increasing 
ls O content on the system is similar to that of increasing Eu content. The band structure calcula- 
tions of the energy gap between t^g and e g bands including the renormalization of this gap due to 
the electron-phonon interaction reveal the physical mechanisms underlying such similarity. 

PACS numbers: 72.80.Ga, 75.30.Wx, 32.10.Bi, 71.27. +a, 71.15.Mb, 75.25.Dk 



I. INTRODUCTION 

Most magnetic oxides are characterized by a strong 
interplay of electron, lattice, and spin degrees of free- 
dom giving rise to multiple phase transitions and differ- 
ent types of ordering. The phase transitions are often ac- 
companied by the formation of different inhomogeneous 
states. In such situation, the oxygen isotope substitution 
provides a unique tool for investigating inhomogeneous 
states in magnetic oxides, which allows studying the evo- 
lution of their properties in a wide range of the phase 
diagram. Sometimes, particularly if a system is close to 
the crossover between different states (usually leading to 
a phase separation), the isotope substitution can lead to 
significant changes in the ground state of the system.—. 

A good example of such phenomena is provided by 
cobaltites. These perovskite cobalt oxides have attracted 
a special interest owing to the possibility of the spin-state 
transitions (SST) for the Co ions induced by tempera- 
ture or doping 2 ^ 8 - and the related phase separation phe- 
nomenal— The effect of 16 O — > 18 isotope substitution 
on the properties of (Pri_. y Euj / )o.7Cao.3Co03 cobaltites 
(0.12 < y < 0.26) was studied earlier in our paper Ref.ll7l. 
It was found that with increasing Eu content, the ground 
state of the compound changes from a "nearly-metallic" 
ferromagnet (ferromagnetic metallic clusters embedded 
into an insulating host) to a "weakly- magnetic insula- 
tor" at y < y cr ps 0.18, regardless the isotope content. A 
pronounced SST was observed in the insulating phase (in 
the samples with y > y cr ), whereas in the metallic phase 



(at y < y cr ), the magnetic properties were quite differ- 
ent, without any indications of a temperature-induced 
SST. Using the magnetic, electrical, and thermal data, 
we constructed the phase diagram for this material. The 
characteristic feature of this phase diagram is a broad 
crossover range near y cr corresponding to a competition 
of the phases mentioned above. The 16 O —> 18 substitu- 
tion gives rise to an increase in temperature Tss of the 
SST and to a slight decrease in the ferromagnetic (FM) 
transition temperature Tfm- 

However, a number of problems important for under- 
standing the physics of the systems with spin-state tran- 
sitions have not been touched upon in the study reported 
in Ref. Il~7l . The most important question is the relation 
between the changes caused by varying the composition 
(increase of concentration y of the smaller rare-earth ions 
Eu) and that due to isotope substitution, and the phys- 
ical mechanism underlying these changes. We see, c/. 
phase diagram in Ref. [171 and in Fig. [1~2"1 below, that there 
exists some correlation between these changes, but the 
situation is not so simple: thus, in the right part of the 
phase diagram the spin state transition increases both 
with the increase of Eu content y and with the increase 
of isotope mass (going from 16 to 18 0). At the same 
time, in the left part of this phase diagram the effect of 
the increase of Eu content and of increasing oxygen mass 
on the phase transition (in this case transition to a nearly 
ferromagnetic state) is just the opposite: an increase in 
Eu content leads to the decrease of Tfm, but the increase 
of oxygen mass - to the increase of it. 
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Another important open question concerns the behav- 
ior of separate phases in the regime of phase separation. 
There are many different correlated systems, in which 
phase separation was detected in some range of compo- 
sitions, temperatures, external fields, etc. Typically, the 
measured transition temperatures in this case changes 
e.g. with doping. However, it often remains unclear 
whether this change is the effect occurring in separate 
regions of different phases, or is just the result of averag- 
ing out over the inhomogeneous system. To answer these 
questions, we now carried out the detailed study of the 
behavior of the (PrEu)CoC>3, using as a tool the possibil- 
ity of fine tuning the properties of the system by partial 
isotope substitution. This partial substitution plays in ef- 
fect the same role as doping, external pressure, etc. The 
obtained results established the possibility of "rcscaling" 
the changes in the system with doping and with isotope 
substitution and allowed us to clarify the questions for- 
mulated above. 

As to the second question formulated above, just the 
possibility of fine tuning the properties of the system in- 
side the region of phase separation, provided by partial 
isotope substitution, allowed us to study separately the 
behavior of different phases within this phase separated 
regime - the possibility, which would be very difficult to 
get by other means. Our results obtained in this way 
demonstrate that not only the average critical tempera- 
tures change with doping and with isotope substitution, 
but also "individual" transition temperatures (tempera- 
ture of ferromagnetic transition in more metallic regions 
and the temperature of the spin-state transition in more 
insulating parts of the sample) do change with chemical 
and isotope composition. 

As to the main, the first question formulated above, 
about the mechanisms governing the change of proper- 
ties of the system with chemical and isotope composition, 
the experimental findings reported in the present paper 
provided us an opportunity to formulate a realistic the- 
oretical model clarifying the mechanisms underlying the 
pronounced isotope effects in cobaltitcs exhibiting spin- 
state transitions. The theoretical analysis demonstrates 
that the main factor is the change of the effective band- 
width with the change both of chemical and isotope com- 
position. The opposite trends in two parts of phase di- 
agram mentioned above find natural explanation in this 
picture. 

To analyze the effects of the partial oxygen isotope 
substitution for the doped cobaltites in the crossover 
region of the phase diagram, we have prepared a se- 
ries of oxide materials with nearly continuous tuning 
of their characteristics. This allows us to trace the 
evolution of relative content of different phases as a 
function of 18 0/ 16 ratio. Note that there were only 
few investigations of this kind, one of which we under- 
took earlier for (Lai_ 2/ Pr 3/ )o.7Cao.3Mn03 manganitesji& 
Here, the pronounced isotope effect manifesting itself 
in (Pri_. y Euj / )o.7Cao.3Co03 cobaltites provides indeed 
a unique possibility to address the problems discussed 



above through the use of the partial oxygen isotope sub- 
stitution. 



II. SAMPLES 

Polycrystalline (Pri_j / Eu y )o.7Ca .3Co03 samples were 
prepared by the chemical homogenization ("paper syn- 
thesis" ) method^ through the use of the following opera- 
tions. At first, non-concentrated water solutions of metal 
nitrates Pr(N0 3 ) 3 , Eu(N0 3 ) 3 , Ca(N0 3 ) 2 , and Co(N0 3 ) 2 
of 99.95% purity were prepared. The exact concentra- 
tion of dissolved chemicals was established by gravimetric 
titration and, in the case of Co-based solution, by means 
of potentiometric titration. The weighted amounts of 
metal nitrate solutions were mixed in stoichiometric ratio 
and the calculated mixture of nitrates were dropped onto 
the ash-free paper filters. The filters were dried out at 
about 80° C and the procedure of the solutions dropping 
was performed repeatedly. Then, the filters were burned 
out and the remaining ash was thoroughly ground. It 
was annealed at 800°C for 2 hours to remove the car- 
bon. The powder obtained was pressed into the pellets 
and sintered at 1000°C in the oxygen atmosphere for 100 
hours. Finally, the samples were slowly cooled down to 
room temperature by switching off the furnace. 

Samples were analyzed at room temperature by the 
powder X-ray diffraction using Cu K a radiation. All 
detectable peaks were indexed by the space group, 
Pnma. According to the X-ray diffraction patterns, 
all (Pri_.yEu w )o.7Cao.3Co03 samples were obtained as 
single-phase crystalline materials. 

We prepared a series of ceramic cobaltite samples with 
the degrees of enrichment by 18 equal to 90%, 67%, 
43%, 17%, and 0%. These values were determined by 
the changes in the sample mass in the course of the iso- 
tope exchange and by the mass spectrometry of the resid- 
ual gas in the oxygen exchange contour. The samples 
were annealed in the appropriate 16 0- 18 gas mixture 
at 950 °C during 48 h at total pressure of 1 bar. The 
Eu doping of the samples was chosen to be near and on 
the both sides of the crossover doping level y cr : a com- 
position with the high Eu content exhibiting the SST, 
a low-Eu composition corresponding to the nearly ferro- 
magnetic state, and the sample at the phase crossover, 
where both Tss and Tfm, were observed. For all sam- 
ples, we measured the temperature dependence of the 
real x'(^) an d imaginary \'' '(T) parts of the ac magnetic 
susceptibility and electrical resistivity p(T). The resis- 
tivity of the samples was measured by the conventional 
four-probe technique in the temperature range from 4.2 
to 300 K . The measurements of ac magnetic suscep- 
tibility x(T) were performed in ac magnetic field with 
a frequency 667 Hz and an amplitude of about 0.1 Oc. 
Based on these measurements, we were able to analyze 
the evolution the sample properties with the change Eu 
and 18 content. 
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Figure 1: (Color online) Temperature dependence of magnetic 
susceptibility for (Pri_ H Euj,)o.7Cao.3CoC>3 with y = 0.2. 



III. EXPERIMENTAL RESULTS 

1. For the samples with y > y cr (Eu content y = 0.20), 
the material correspond to a "weakly-magnetic insula- 
tor", in notation of Ref. [l?]. The x'(T) curves give clear 
indications of a spin-state transition (SST) at T$s man- 
ifesting itself as a peak in the x'(T), see Fig. Q] The 
value of T$s increases with the 18 content. The high- 
temperature phase is a paramagnet and a relatively good 
conductor, see Fig.[2j The low-spin (LS) insulating phase 
is dominant below the crossover to a LS state occurring 
at about 100 K. The increase of \' a t low temperatures 
is most probable caused by an incomplete transition, af- 
ter which there may remain small magnetic (and pre- 
sumably more conducting) clusters immersed into the 
LS state bulk insulator. As a result of this crossover 
to a LS state, the electrical resistance R increases by 
10-12 orders of magnitude, which can be treated as a 
metal-insulator transition (MI), see Fig. (2). In these 
compounds, the metal-insulator transition is accompa- 
nied (or caused) by the SST. The ls O content does not 
produce a significant effect on R(T) at low T, although R 
is larger in the samples with heavy oxygen. The value of 
the metal-insulator transition temperature Tmi can be 
determined from the logarithmic derivative of R{T), see 
inset of Fig. [2] The mass dependence of Tmi correlates 
well with the isotopic shift of T$s, see Fig. [3J Accord- 
ing to the calculation of isotopic constant T ~ Af~ Q ; 
a = -dhxT/dhiM = — (AT /AM) (M/T) , we have the 
value ass and olmi = —(1-5 ± 0.07). The increasing in 
the oxygen mass promotes the development of the LS 
state. 

2. The most important results are obtained for the 
samples with Eu content y ~ y cr (y = 0.14 and 0.16). 
They correspond to a wide concentration range of the 
phase separation. 

For the sample with y = 0.14, we observe in x'(T) 
curves a feature corresponding to a steep increase of the 



Figure 2: (Color online) Temperature dependence of electrical 
resistivity for (Pri_ y Eu !/ )o.7Cao.3Co03 with y = 0.2. 
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Figure 3: (Color online) Isotopic shift of the characteristic 
temperatures of the spin-state Tss and metal-insulator Tmi 
transitions at the ie O— > ls O substitution for the sample with 
2/ = 0.2. 

magnetization on cooling at 60-70 K, see Fig. 2] This 
behavior is caused by the FM phase arising in these sam- 
ples. In addition, in the temperature dependence of resis- 
tivity, we see a steep increase of resistivity characteristic 
of the metal-insulator transition similar to that observed 
for the samples with y = 0.2, see Fig. [5] The transition 
temperature corresponds to the minimum in the logarith- 
mic derivative of the resistivity, see inset of Fig. [5J This 
means that here we deal also with the change in the rela- 
tive content of metallic and insulating phases suggesting 
the existence of the regions corresponding to the low-spin 
insulating state and the correlation between the metal- 
insulator and SS transitions. At the same time we do not 
observe any clear indications of the SST in the temper- 
ature dependence of the magnetic susceptibility. Thus, 
the samples with y = 0.14 turn out to be at the bound- 
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Figure 6: (Color online) Isotopic shift of the characteristic 
Figure 4: (Color online)Temperature dependence of magnetic temperatures of the spin-state T ss , metal-insulator T M i and 
susceptibility for (Pr 1 _ v E % ) .7Cao. 3 CoO 3 with y = 0.14. ferromagnetic T FM transitions for the sample with y = 0.14. 




Figure 5: (Color online) Temperature dependence of electrical 
resistivity for (Pri_ y Eu H )o.7Cao.3Co03 with y = 0.14 (solid 
line) and y = 0.1 with ie O (dashed line) and 18 (dash-dot 
line). 



Figure 7: (Color online) Temperature dependence of magnetic 
susceptibility for (Pri_j / Eu y )o.7Cao.3Co03 with y = 0.16. 



ary of the phase separation range. The increasing in the 
oxygen mass favors the LS state as well as suppression of 
the FM phase and of the metallization. 

The values of the isotopic constant calculated for the 
SS and FM transitions in this sample are ass. mi = 
-(1.7 ± 0.06) and a FM = (0.5 ± 0.05), respectively. For 
this sample, the isotope shifts of the characteristic tran- 
sition temperatures arc illustrated in Fig. [6l 

For the samples with y — 0.16, the effect of the partial 
oxygen isotope substitution by 18 manifests itself even 
more clearly. Here, we observed the features characteris- 
tic both of FM and SS transitions. On the one hand, the 
temperatures dependence of %' exhibits a steep growth at 
60-70 K similar to that in the sample 0.14 indicating the 
existence of the FM transition. On the other hand, in of 
x'(T) , we observed a clearly pronounced peak, which can 
be attributed to SST at T ss (Fig. [7]). x'(T) curves also 



demonstrate that the transition to the LS state gradu- 
ally disappears with the decrease of the oxygen mass. In 
samples with small 18 content (< 17%), this transition 
is hardly seen due to the gradual transformation from 
the LS to FM state when the oxygen mass decreases. 
Note also that for y = 0.16 (as well as in the samples 
with y = 0.14) the temperature Tfm decreases with the 
growth of the average mass of oxygen. 

In addition, the Curie- Weiss temperature 9 in the for- 
mulas for the inverse magnetic susceptibility x _1 (T) con- 
siderably decreases with the increase of the average mass 
of oxygen isotopes (see Fig. [5]). For the samples with the 
largest oxygen mass, we have 6 = 184 K, whereas for the 
samples with the 16 0, 9 = 45 K. This phenomenon may 
be related to the transition from the antiferromagnetic 
interaction to ferromagnetic one. 

For sample with y = 0.16 as well as in the samples with 
y = 0.14 and 0.2, we have found the growth of resistivity 
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Figure 8: (Color online) Temperature dependence of inverse 
magnetic susceptibility for (Pri_j,Eu K )o.7Cao.3Co03 with y = 
0.16. 
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Figure 9: (Color online) Temperature dependence of electrical 
resistivity for (Pri_j / Eu H )o.7Cao.3Co03 with y = 0.16. 

at decreasing temperatures (by 5-7 orders of magnitude) 
and the metal-insulator transition in the vicinity of 70 K 
(Fig. [5]). Both the resistivity and magnetic susceptibil- 
ity data give clear indication that this sample is in the 
phase-separation range. The temperatures of MI and SS 
transitions increase with the average mass oxygen (see 
inset of Fig. In the sample with y = 0.16, we see the 
same general tendency, namely the growth of T$s and the 
decrease of Tfm with the growth of the average oxygen 
isotope mass. Here, the values of the isotope constant are 
&SS,mi = — (I-? ^ 0.06) and apM = (0.4 ± 0.1), respec- 
tively; the isotope shifts of the characteristic transition 
temperatures are illustrated in Fig. [TO] 

3. Finally, the samples with y < y cr (Eu content 
y = 0.10) fall into the range of "nearly-metallic" fer- 
romagnet. In Ref. [l?], it was shown that at T < Tfm the 
compositions with a low Eu content correspond to the 
domains of the metallic ferromagnetic phase embedded 
in a weakly magnetic nonconducting matrix. 



Figure 10: (Color online) Isotopic shift of the characteristic 
temperatures of the spin-state Tss, metal-insulator Tmi and 
ferromagnetic Tfm transitions for the sample with y = 0.16. 
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Figure 11: (Color online) Temperature dependence of mag- 
netic susceptibility for (Pri_ a Euj,)o.7Cao.3Co03 with y = 0.1. 



According to the temperature dependence x'(T) plot- 
ted in Fig. [TT] the magnetization steeply increases on 
cooling at about 60 — 70 K (i.e. the ferromagnetic phase 
arises). With the growth of the average oxygen mass, the 
value of x'(T) decreases at low temperatures. Here, we 
have Tfm{ 1s O)< Tfm( 16 0) and the maximum isotope 
shift of Tfm does not exceed 2-3 K. 

The electrical resistivity for this samples (y = 0.10) 
increases at low temperatures, even when the ferromag- 
netic phase arises. This behavior is quite similar to that 
observed in the samples with y = 0.14 and y = 0.16. 
Note that the electrical resistance in the samples with 
18 is higher than in those with 16 0. However, with the 
decrease of Eu content down to y = 0.10, the MI pe- 
culiarity in the resistance is suppressed. Such behavior 
R(T) for samples with y = 0.10 can be compared with 
that for y = 0.14 samples, see Fig. and inset of this 
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figure. The resistivity for samples y = 0.10 corresponds 
to a more smooth curve than for samples with y = 0.14, 
although the regular course of the temperature depen- 
dence remains nearly unchanged. Thus, the samples with 
y = 0.10 do not become really metallic but their behavior 
differs from the behavior of the samples with y = 0.14 
(they do not exhibit indications of SST). Therefore, we 
argue that in the phase diagram, the composition with 
y = 0.10 lies outside the crossover region, on the left- 
hand side of it. 



IV. DISCUSSION OF EXPERIMENTAL 
RESULTS 

The obtained data for the temperatures of the phase 
transitions can be presented in the form of the phase di- 
agram (Fig. [T2j) . It illustrates that with the growth of 
Eu content the system transforms from the nearly fer- 
romagnetic metal to the LS insulating state undergoing 
LS — > IS spin-state transitions . Between these states, we 
have a broad crossover region corresponding to the phase 
separation. Indeed, the simultaneous observation both of 
Tfm and T$s in the samples with Eu content y = 0.14 
and 0.16 is a clear evidence of the phase separation in 
the system. 

These samples also provide a spectacular illustration of 
the effect related to the variable content of 18 0. In par- 
ticular, for the samples with y = 0.16, the temperature 
dependence of magnetic susceptibility (Fig. [3J exhibits at 
high values of 18 O content a pronounced feature corre- 
sponding to the SST. With the lowering of 18 content, 
this feature becomes weaker and disappears below 17% 
of 18 0. Hence, we see that the change on the average 
oxygen mass can drastically affect the phase composition 
of the cobaltite samples. 

The oxygen isotope substitution 16 — > 18 shifts the 
phase equilibrium toward the insulating state. For the 
heavier isotope, the spin-state transition temperature 
Tss grows, while the ferromagnetic transition tempera- 
ture Tpm decreases. The change of average oxygen mass 
is a unique tool for investigating special properties of 
phase separation in cobaltites near the crossover between 
the FM and LS phases. We also see that the effect of in- 
creasing 18 content on the system is similar to that of 
increasing Eu content. 

The analysis of these results for different chemical and 
isotope composition demonstrates that the effect of the 
increase in Eu content y and of average oxygen mass are 
qualitatively similar. One can rescale the dependence of 
the transition temperature on both parameters using the 
combined variable y + 0.01a;, where x is the relative con- 
tent of 18 0, as shown in Fig. [13] Both the temperatures 
of SS and FM transitions depend almost linearly on this 
combined variable. We see that the change of Eu content 
by 1% is equivalent to the change of the isotope content 
by 100 % . This actually the most important result of the 
present study. The theoretical analysis of these results is 
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Figure 12: (Color online) Phase diagram of 

(Pri- y Eu !/ )o.7Cao.3Co03 compound with ie O and 90% 
18 0. PM, FM, IS, and LS stand for the paramagnetic, ferro- 
magnetic, intermediate-spin, and low-spin state, respectively. 
The hatched area corresponds to the phase-separated state 
(PS). 



given in the following sections. 



CALCULATION DETAILS 



To explain the composition and isotope dependence of 
the properties of our system (see Fig. [TJ]), especially the 
similar dependence of the SST temperature and the tem- 
perature of the FM transition illustrated in Fig. [TjJJ we 

based predom- 



propose a realistic model (Section IVIIJ) 
inantly on the change of the electron bandwidth with 
chemical and isotope composition. Some input, as well 
as the estimates of relevant parameters are taken from 
the ab initio band structure calculations for the limiting 
"pure" compositions corresponding to y = (PrCoOs) 
and y = 1 (EuCo0 3 ). 

The crystal structure of PrCoOa obtained in Ref. 
for T = 300 K was utilized in those calculations. For 
EUC0O3, the lattice parameters were taken from Ref. I23I . 
The exact atomic positions for EUC0O3 arc unknown, 
therefore we used the same positions as for PrCoOa (with 
the correct unit cell volume for E11C0O3). The splitting 
between different one-electron energy levels Acf was cal- 
culated within the local density approximation (LDA) in 
the framework of the linear muffin-tin orbitals method 
(LMTO)i24 Partially filled, but physically unimportant 
4/ states of the Eu and Pr were treated as frozen^ 

The Brillouin-zone (BZ) integration in the course of 
self-consistency iterations was performed over a mesh of 
144 k-points in the irreducible part of the BZ. 
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Figure 13: (Color online) Plots illustrating the combined ef- 
fect of the Eu and 18 O doping on the temperature of spin-state 
transition Tss and on the ferromagnetic transition tempera- 
ture Tfm- Note the difference of the temperature scales in 
panels (a) and (b). 



VI. DOPING DEPENDENCE: LDA RESULTS 

There arc different ways to estimate the temperature of 
the spin-state transition with the use of the band struc- 
ture calculation. The most direct way is to calculate 
the total energies of different spin statesi 20 i 21 . However, 
in the case of the doped system this would require very 
large supercells. Moreover, currently,we have no single 
commonly accepted model, which can explain all exper- 
imental facts. Various combinations of a static or dy- 
namic order of the different spin states are discussed in 
the literaturej 25 ' 39 " — This is the reason why we chose an 
alternative approach. 

The energy of any of the spin states depends on two im- 
portant parameters: the single-electron energy difference 
Acfs between the highest ti g and the lowest e g levels, 
and the intra-atomic Hund's rule exchange coupling Jh ■ 
The Hund's rule energy Jh is an atomic characteristic 
and does not change appreciably cither with the Eu dop- 



ing or with the isotope substitution. So, to investigate 
the dependence of the spin-state transition temperature 
on the doping or isotope substitution in the first approx- 
imation one can focus on the study of the single param- 
eter, A C fs- If A C fs is large enough (A C fs > 2Jjj), it 
will be energetically favorable to localize all electrons in 
the low- lying t 2g subshcll of Co 3+ , i.e. the system will 
be in the low-spin state. With decrease of the crystal- 
field splitting, some electrons can be transferred to e g 
subshcll, which allows the system to gain the exchange 
energy, since there will be more electrons with the same 
spin. As we will show below both the isotope substitu- 
tion and the doping can be related to the crystal-field 
splitting. 

In order to estimate Acfs, we used the Wannier func- 
tion projection procedure proposed in Ref. I2H which al- 
lows us to project the full-orbital band Hamiltonian onto 
the subspace of few states (five d states of Co). With the 
Fourier transformation one obtains the Hamiltonian in 
the real space, from which the splitting between highest 
in energy t 2g and lowest e g can be easily calculated. For 
PrCo0 3 , we obtain A C fs = 2.07 eV. 

The total and partial densities of states (DOS) ob- 
tained for PrCoC>3 in the LDA calculations are presented 
in Fig. [L4j In the octahedral symmetry, the 3c? states 
of Co arc split into t 2g and e g subbands. In the LDA, 
the valence band is mostly formed by the Co— t 2g states, 
while the conduction band is determined by the Co— e g 
states. The O— 2p band is located in the energy range 
from -7 to -1.5 eV. 

The DOS for EUC0O3 is qualitatively very similar and 
corresponding calculations of Acfs result in the value of 
2.14 eV. The increase of the t2 g — e g excitation energy on 
going from PrCoOs to EUC0O3 is caused by two factors. 
The first is the lanthanide contraction: the substitution 
of the large Pr 3+ by smaller Eu 3+ ions leads to some 
decrease of the Co-0 distance and to the corresponding 
increase in the p — d hybridization, which leads to an in- 
crease of the difference between the centers of the t2 g and 
e g bands. The second effect is related to the decrease of 
the effective widths of t 2g and e g energy bands with the 
corresponding increase of the energy gap between them. 
This narrowing of energy bands on going from PrCoOa to 
EUC0O3 is also related, in effect, with the lanthanide con- 
traction, because of which the tilting of C0O6 octahedra 
increases and the Co-O-Co angle and the corresponding 
bandwidth decrease on going from PrCoOs to EUC0O3. 
Both these effects finally lead to the increase of Acfs 
with Eu content, which leads to the enhanced stabiliza- 
tion of of the LS state of Co 3+ (see the more detailed 
discussion of these effects in Section |VII[) . 

This change of the crystal-field splitting (CFS) results 
in the modification of the spin-state transition temper- 
ature, since this transition is due to the competition of 
the Hund's rule exchange coupling Jh and the CFSi2£ 

It was found in Refs. [T3 and (see also 

Fig. [T2|) that the change of the Eu content y in 
(Pri_j / Eu y )o.7Cao. 3 Co03 by 0.02 leads to the change of 
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Figure 14: (Color online) The total and partial densities of 
states (DOS) for PrCoC>3. The Fermi energy corresponds to 



the spin-state transition temperature by about 14 K. 
In the first approximation, it is possible to neglect the 
presence of Ca and interpolate the change of the CFS 
for the complex system like (Pri_j,Eu y )o.7Cao.3Co03 us- 
ing the values of the CFS for y = 1 (PrCoOs) and 
y = (E11C0O3). Such a linear interpolation predicts the 
change of the spin-state transition temperature by 16 K, 
if y changes by 0.02, which is in an excellent agreement 
with the exp er iment 1 1 7 1 2 7 



VII. ISOTOPE SUBSTITUTION: MODEL 
RESULTS 

The isotope substitution does not change the chemical 
properties of the ions such as the oxidation numbers or 
bonding energies. However, it affects the crystal lattice 
through the modification of the phonon spectra. Below, 
following the approach of Ref. [28|, we demonstrate the 
effect of this modification on the electronic and magnetic 
properties and hence on the spin-state transition. 



In the tight-binding model, the band spectra of a solid 
is determined by the on-site ionic energy levels ef lm 
and the hopping matrix elements between different sites 
^ ,mm ^ -phe ionic energies e™ /m obviously do not de- 
pend on the mass of the ions, being determined by the 
quantum numbers and the intra-atomic Coulomb and ex- 
change interactions. The hopping parameters depend on 
the type of the orbitals (s,p, d, /), bonding type (jr, a, 6), 
and the distance between ions, u. According to the fa- 
mous Harrison parametrizatio n 29 i 30 in the absence lattice 
vibrations the hopping integrals between e.g. p orbitals 
of the oxygen and transition metal <i-orbitals equal to 



t 



C 



pd 



pdm 



(l) 



where coefficients C p d m depend on the type of the bond- 
ing and can be different for different metals and lig- 
ands^^ 

The static version of Eq. (Q} can be generalized tak- 
ing into account the presence of lattice vibrations, i.e. 
phonons, which depend on the ion masses. The mean pd 
hopping matrix element can be calculated as 



(t 



1 



pd) 



u +v 



pdm 



2« Ju -v 



du 



a 



pdm 



1 



1 



6v V (uq — v) 3 (uq + v) ? - 



(2) 



where v = y (5u 2 ) is the mean square displacement 
from the equilibrium position uo due to phonons. Since 
v/uq <C 1 one may simplify last equation expanding it to 
series to the 4th order 



(t P d) 



pdm 



1 



3 \uq 



o((f) 1 ) 



(3) 



In the static limit v — > 0, the last formula coincides 
with Eq. (P). 

In the Debye model at zero temperature, the mean 
square displacement is written as^ 



m 2 



1 



4fc_B^D m 



(4) 



Here m is the mass of vibrating ions and 9d is the Debye 
temperature. Due to different mass the mean square dis- 
placement in the compounds enriched by 16 O or 18 O will 
be different. The Debye temperature for the very similar 
system, LaCo03 was found to be ~600 K<2£. One may 
use this value to estimate v in (Pri_ l( Euj / )o.7Cao.3Co03. 
Then the mean square displacement for 18 is vig = 
0.100 A and v w = 0.107 1 for 16 0. According to Eq. ©, 
this leads to the decrease of the effective bandwidth in 
going from 16 to 18 0. Qualitative explanation of this ef- 
fect is presented in Fig. 115b . For strong electron-phonon 
coupling the same effect could be attributed to polaron 
band narrowing depending on the isotope mass. 
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The temperature of the spin-state transition depends 
on the energy difference between t 2g and e g subbands, 
which is defined by the widths of corresponding bands 
and the positions of their centers. We start with the 
study of the bandwidths dependence on the ligand ion 
mass, see schematic illustration in Fig. [T5b . 

In order to calculate the change in the bandwidth 
caused by the 16 — > 18 substitution, one needs to know 
the hopping integrals, which depend on two unknown 
parameters C p d m and uq. The C p d m coefficients can be 
in principle evaluated as it is prescribed for instance in 
Ref. [29|. However, for the better precision we calculated 
Cpdm parameters from the LDA t 2g and e g bandwidths in 
pure PrCoC>3. The equilibrium Co-0 distance uq in its 
turn can be evaluated from the actual crystal structure 
data for EUC0O3 and PrCoOs. 

Finally, performing all these calculations, one gets that 
the t2 g bandwidth decreases on 22 K when one substi- 
tutes 16 O by 18 O. The decrease of the e g bandwidth is 
two times larger and equals 44 K. In effect, the minimum 
energy of the t 2g — > e g transition increases by 33 K. Hence 
the spin-state transition temperature due to the change 
of the bandwidth must increase by the same value go- 
ing from 16 O to 18 O in qualitative accordance with the 
experiment 

Let us consider the second mechanism, which affects 
the spin-state transition and which is related to the de- 
pendence of changes in the centers of gravity of cor- 
responding bands (i.e. crystal-field splitting) with the 
isotope substitution. It turns out the this effect coun- 
teracts the first onc(changc of the effective t 2g and e g 
bandwidths), but numerically this second effect is much 
smaller, see below. Generally speaking, there are two 
main contributions to the crystal-field splitting, Aqfs- 
One comes from the Coulomb interaction of the 3d elec- 
trons with the negatively charged ligands, another is due 
to the hybridization between d orbitals of metal ions 
metal and p orbitals of the ligands . 1 35 For most ox- 
ides of 3d transition metals both terms are acting in "the 
same direction" , resulting to the same sequence of lev- 
els^ That is why so crude approaches as the atomic 
sphere approximation (ASA)2£ often used in the ab initio 
calculations provide quite precise band structure in most 
cases. The effect of the Coulomb term can be omitted or 
effectively incorporated into the kinetic energy contribu- 
tion. Below we will follow the same strategy by consid- 
ering the kinetic energy only, keeping in the mind that 
the Coulomb contribution can be taken into account via 
the parameters renormalization. 

In the second order of the perturbation theory, the 
crystal-field splitting between t 2g and e g subbands is 
written as 

t 2 - t 2 

&CFS - T , (5) 

Act 

where Act is the charge-transfer energy (which corre- 
sponds to the d n p 6 — > d n+1 p 5 transition), t p da and t p d-K 
are the hopping matrix elements for different types of the 



bonds. 

Since the average hopping {t p d} decreases on going 
from 16 to 18 according to Eqs. © and the 
crystal-field splitting should also decrease as follows from 
Eq. (|SJ). As a result, this contribution should lead to the 
opposite tendency: decrease of the spin-state transition 
temperature going from 16 to 18 0. However, this effect 
does not exceed few kelvins, at the realistic values of the 
charge-transfer energy in cobaltites^ 

Note also that we estimated here the changes in the 
distance between the edges of t 2g and e g subbands. How- 
ever, at finite temperatures, we should have not only the 
transitions between the band edge but just from one sub- 
band to another. Such a temperature-induced smearing 
could diminish somehow our estimates of the isotope ef- 
fect in SST. 

In Fig. [T3"b . we see that the isotope effect for Tfm, 
being much weaker (than for T55), is of the opposite 
sign. Nevertheless, there is the same similarity between 
the effects of the Eu content and the oxygen isotope sub- 
stitution. This is in agreement with our expectations, 
because the ferromagnetism of the low-Eu doped sam- 
ples with metallic clusters should be stabilized by the 
double-exchange mechanism, according to which Tfm is 
proportional to the effective bandwidth of the itinerant 
electrons. This bandwidth decreases for the heavier iso- 
tope, that is why Tfm goes down at the 16 — > 18 sub- 
stitution as well as at increasing Eu content, see Fig. 113b. 
A schematic illustration of the mechanism underlying the 
oxygen isotope effect discussed above is given in Fig. [15] 



VIII. CONCLUSIONS 

Experimental studies carried out for 
(Pri_ y Eu y )o.7Cao.3Co03 cobaltites with varying 
isotope substitution of 16 by 18 demonstrated that 
there exists a strong similarity in the changes caused 
by the chemical composition (increasing Eu content) 
and those arising from the oxygen isotope substitution. 
The chemical composition y ~ 0.1 — 0.2 was chosen 
because in this range a crossover occurs between the 
ferromagnetic near-metallic state with magnetic Co ions 
to the nonmagnetic insulator with the low-spin Co 3+ 

(tg fl e2 J S = 0),seeFig.H2 

The main experimental conclusion presented in Fig. [T3] 
is that one can rescale the behavior of this system. The 
dependence of the spin-state transition temperature and 
of the ferromagnetic transition temperature on content 
y of Eu and on the content x of the heavier isotope 18 O 
can be represented by the same almost linear plot as func- 
tion of the combined variable y + O.Olx. This means, for 
example, that the change of Eu content y by 0.005 is 
equivalent to the substitution of 50% of 16 by 18 0. In 
addition, it clearly demonstrates that not only the aver- 
age transition temperatures change with doping and with 
isotope substitution, but also the transition temperatures 
for each separate phase vary with chemical and isotope 
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relative content of O 

(b) 

Figure 15: (Color online) Schematic illustration of the ef- 
fects of oxygen isotope substitution (a) on the effective hop- 
ping integral t e ff and (b) on the effective gap A e // be- 
tween the e g and ti g bands and on the effective width W e ff 
of the e 9 band. Solid curve in panel (a) depicts the hop- 
ping integral t as function of interion distance u accord- 
ing to Eq. JTJ. Effective hopping integrals t e ff are deter- 
mined by averaging t(u) over the interval given by the mean 
square displacement v of ions due to the lattice vibrations 
Q (crudely, t e ff = l/2[t(«o + v) + t(uo — v)] for each iso- 
tope). We see that i e //( 16 0) > i e //( 18 0). As a result, 
A e// ( ls O) > A e// ( 16 0) and hence T ss ( 18 0) > T ss ( 16 0). At 
the same time, W eff { 18 0) < W eff ( 1(i O) and hence T FM ( 18 0) 
< TW( 18 0). 



composition. 

Based on this similarity between the role of chemical 
and isotope composition for the spin-state transition and 
for the transition to the ferromagnetic state at a smaller 
Eu content, we propose a theoretical explanation of the 
isotope effect in these transitions. We investigate the cor- 
responding changes and estimate the relevant parameters 
using the ab initio band structure calculations. These 
results together with the analytical model allow us to ex- 
plain the observed behavior. In particular, the isotope 
effect both in the spin-state and ferromagnetic transi- 
tions is interpreted in terms of the change in the cor- 
responding widths of the d bands occurring due to the 
electron-phonon renormalization, which depends on the 
atomic masses of the respective isotopes. 

All this demonstrates once again that the oxygen iso- 
tope substitution is a powerful tool for revealing salient 
features in the behavior of strongly correlated magnetic 
oxides. 

Summarizing, we can say that using this approach, we 
established, first, that in the aforementioned crossover 
range, we can clearly distinguish two coexisting phases, 
nearly insulating exhibiting a spin-state transition and 
nearly metallic "ferromagnetic" , with the different be- 
havior of the transition temperatures. Second, we have 
found that these transition temperatures depend almost 
linearly on the content of the heavy oxygen isotope, which 
is a non-trivial observation clearly demonstrating that 
the electronic structure could be effectively controlled by 
isotopes. Third, based on these observations and using 
the parameters deduced from our band-structure calcu- 
lations, we put forward a simplified model capturing the 
main physics of the isotope effect in the systems with 
spin-state transitions and quantitatively describing the 
experimental data. 
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